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FOREWORD 

The  research  described  in  this  report  was  performed  during  the  1993  fiscal  year  as  part 
cf  an  effort  to  improve  radar  classification  ana  identification  capabilities  for  non- 
cooperative  targets  This  problem  continues  to  be  a  primary  goal  or  radar  research 
program'  and  considerable  effort  has  been  expended  within  the  last  few  decades  in 
attempts  to  solve  it  The  accuracy  of  these  methods  depends  upon  avadable  radar 
resolution  and  so,  indirectly,  it  depends  on  radar  bandwidth 

A  novel  technique  for  radar-based  non-cooperative  target  recognition  (NCTR)  lias 
been  created  and  verified  using  a  variety  of  real  and  synthesized  data.  This  method  uses 
ordinary  azimuth-eievation  phase  monopulse  measurements  and  is  appropriate  for 
application  in  existing  (narrow-band)  missile  systems.  An  enhancement  bar  also  oeen 
developed  (using  wtde-band  measurements)  and  offers  significantly  greater  target 
classification  ability  This  synopsis  briefly  describes  this  practical  algorithm  and  lists  its 
required  input  data  as  well  as  displaying  an  example  of  its  ir.ent. 
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INTRODUCTION 


A  problem  of  recurring  interest  to  radar  engineers  concerns  the  classification  and 
identification  of  aircraft  which  do  not  broadcast  their  own  identity  So-called 
non-cooperative  target  recognition  (NCTR)  continues  to  be  a  primary  goal  of  radar 
research  programs  and  considerable  effort  has  been  expended  within  the  last  few  decades 
in  attempts  to  solve  it.  Unfortunately,  to  date,  few  “appropriate”  techniques  (and  even 
fewer  “practicable”  techniques)  have  been  developed. 

The  variety  of  radar-based  methods  that  have  been  suggested  for  NCTR  can  be 
divided  into  two  general  groups'  those  which  make  use  of  some  physical  process  unique 
to  the  target  and  those  which  attempt  to  “image”  the  target  itself.  The  physical  process 
group  includes  methods  based  upon  engine  modulations  (JEM),  nonlinear-contact  effects 
(metal-metal  junctions  on  the  target),  as  well  as  proposed  methods  based  on  skin 
vibrational  modulation  effects  The  imaging  category  encompasses  SAR/ISAR  based 
methods,  down-range  profiling,  high-resolution  angle  imaging,  and  a  variety  of 
techniques  collectively  known  as  “inverse  scattenng  reconstruction”  and  which  actually 
attempt  to  create  only  a  generalized  “image”  of  the  target.  (C  f.,  [References  1  through 
3],  and  references  cited  therein.) 

The  limited  success  of  these  approaches  often  results  from  a  failure  to  properly 
account  for  restrictions  imposed  by  realistic  radar/target  environments.  JEM,  for 
example,  can  be  considered  to  be  the  most  successful  of  all  the  methods  but  may  only  be 
used  when  the  target  orientation  allows  the  engine  turbine  blades  to  be  “seen”  by  the 
radar.  Other  “physical-process”  based  techniques  have  never  been  convincingly  applied 
to  NCTR.  The  shortcomings  in  this  first  category  of  methods  have  motivated  most 
researchers  to  concentrate  their  efforts  on  perfecting  imaging  approaches.  Microwave 
imaging  attempts,  with  radar,  to  do  what  is  usually  accomplished  with  conventional  light 
sources.  However,  since  image  resolution  is  proportional  to  the  aperture-to-wavelength 
ratio,  potential  systems  require  too-large  apertures  for  many  applications  of  interest. 
(Reduction  of  wavelength  cannot  be  done  indiscriminately  since  atmosphenc  attenuation 
and  adverse  weather  conditions  severely  limit  the  usable  range  of  the  imaging  system — 
this  precludes  optical  imaging  systems  from  potential  long-range  NCTR  schemes.)  The 
aperture  can  be  synthetically  increased,  of  course,  but  synthetic  aperture  imaging  requires 
accurate  determination  of  target  aspect  over  time  (an  often  lengthy  and  unreliable  task) 
Full-blown  inverse  scattenng  methods  suffer  even  more  from  restrictions  on  the 
necessary  data  they  not  only  require  that  the  data  be  matched  to  known  target  aspect 
angles  but  they  also  frequently  require  that  these  data  be  collected  over  large  bandwidths. 

Finally,  down-range  target  profiling  has  been  suggested  as  a  classification  alternative. 
(This  has  been  much  ballyhooed  in  connection  with  ultra-wideband  radar  research.)  In  its 
simplest  form,  a  down-range  profile  consists  of  target  scatterer  strength  integrated  over  a 
sequence  of  range-bins  and  mapped  to  the  target’s  down-range  extent  When  target 
orientation  is  known  in  advance,  this  can  yield  crude  target  structure  information  and  may 
be  used  for  rudimentary  target  classification  However,  the  target’s  orientation  is  usually 
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not  known  in  advance  (it  is  usually  estimated  from  a  detailed  target  track)'  and  a 
preferable  scheme  would  yield  target  orientation  as  part  of  the  image  In  addition,  the 
simple  down-range  profile  cannot  really  offer  the  amount  of  information  required  for 
meaningful  target  classification 

In  addition  to  long-range  and  all-weather  requirements  (which  can  be  satisfied  only 
by  radar  systems),  the  chosen  recognition/classificatior.  scheme  must  be  tailored  to  the 
platform  on  which  it  will  be  used  This  consideration  usually  restricts  the  type  of  data 
that  can  be  employed  and  the  computational  engine  that  will  process  these  data  One  of 
the  most  restrictive  environments  are  that  of  the  small,  radar-based  missile  system — and 
it  is  this  problem  that  we  shall  address  below  Typically,  radar  missile  seekers  devote 
most  of  their  resources  to  the  problems  of  guidance  and  control  and  have  little  (if  any) 
ability  to  acquire  data  which  are  not  related  to  these  primary  responsibilities 

Our  present  discussion  demonstrates  how  ordinary  phase  monopulse  tracking  data  can 
be  used  to  create  an  enhanced  down-range  target  profile  that  offers  considerably  greater 
information  than  usual  profiles  Phase  monopulse  data  consist  of  measurements  of  the 
gradient  of  the  phase-front  of  the  field  scattered  from  the  target.  Tacking  information  is 
deduced  from  this  gradient  since,  for  a  point  scatterer,  both  gradient  and  target  bearing  lie 
normal  to  the  spherical  phase-front  [Reference  4],  However,  extended  targets  will  not 
generally  scatter  spherically  symmetric  wave  fronts  and  the  phase  gradient  can  only  be 
expected  to  point  to  the  target  on  average.  We  begin  by  examining  the  statistical 
behavior  of  this  kind  of  data  and  show  how  these  statistics  are  related  to  the  target’s 
cross-range  structure  Then  we  extend  these  results  to  large  bandwidth  signals. 
(WideDand  signals  have  been  proposed  for  use  in  missile  systems  in  an  effort  to  mitigate 
the  so-called  “glint”  problem,  which  we  shall  discuss  in  the  conclusion.)  Finally,  we 
present  some  examples  using  computer  generated  data  that  demonstrate  the  algorithm's 
behavior  under  ideal  conditions. 


PHASE  DERIVATIVE  STATISTICS 


The  original  application  of  phase  gradient  statistics  to  target  structure  estimation  was 
performed  for  narrowband  signals  in  [References  5  and  6].  We  briefly  review  the 
narrowband  results  here  for  completeness  and  because  they  form  the  basis  for  the  related 
wideband  results. 

We  fix  a  coordinate  system  to  the  target  and  consider  a  monostatic  scattenng  situation 
in  which  the  transmitter  and  receiver  are  co-located.  In  the  weak  scatterer  approximation 

the  (scaled)  far-field  response  Fs(x,t),  at  position  x  and  time  r,  to  a  signal  F,(y,t) 
incident  upon  the  target  is  simply  the  sum  of  the  responses  from  individual  elements  of 
the  target  The  scattered  field  can  be  written  as 


Fs(x,r)  =  j  Fr(2y,t-2\x\/c)p(y)dy  ,  (1) 

D 
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where  D  is  the  support  of  the  local  scattcrer  density  function  p( y),.  c  is  the  signal  speed, 
and  the  factor  of  2  accounts  for  the  (two-way)  signal  travel  distance 


In  the  quadrature  model  Fs  is  analytic,  with  real  and  imaginary  parts  forming  a 

Hilbert  transform  pair,  and  we  can  write  Fs(x,t)  =  \Fs(x,t^exp(i<l>{\,t))  where  is 

the  phase  of  the  scattered  field.  Let  n  denote  the  direction  of  signal  propagation  If  u  is 
a  unit  vector  perpendicular  to  n  then  the  directional,  derivative  of  this  phase  can  be  found 
from 


u  V<p(x,r)  =  /m 


F'Jx.t) 


hM 


u  VF,(x,t) 


(2) 


F,(y,t)  must  satisfy  the  wave  equation  and  can  be  wntten  F,(y,t)=F,(n  y  ct). 
Substitution  yields 


f 

u  -W(p(\,t)  =  1m 

V 


|Fs(x,f)li 


J 


^rrP(y)/r/'(2W+2"  y ~efKv 

N  ) 


(3) 


where  F',(^)  =  dF,(^)/ci^ 

Since  thi  coordinate  system  is  fixed  to  the  target,  a  measured  value 
gt  =  g(n,)  =  u  V(p(x„r)  will  depend  upon  the  target  orientation  through  the  “weighting” 
function  p(y).  Moreover,  in  typical  radar  situations  the  signal  has  high  frequency 

content  (-10  GHr)  and  gt  vanes  rapidly  with  only  minute  changes  in  target  aspect  (It  is 
known,  for  example  that  measurements  from  targets  flying  “straight  and  level  will 

exhibit  significant  /anation.)  For  this  reason,  we  consider  p(y)  to  be  a  random  process 

and  examine  the  statistical  behavior  of  g,. 
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The  narrowband  case  When  we  set  F,(y,t)  =  exp[ik[n  y-cr))  where  k  =  (olc  is  the 
wave  number  associated  with  the  time  harmonic  incident  field  of  frequency  co ,  we  obtain 


u  V<p(x,t)  =  Im\ 


Fi(x,t)e‘ 


If — J  -^p(y)«p(i2A"  y)dy 


|fs(*.0|  D  lxl 


(4) 


Define  G(x,r)  =  [G(x,r)|e,'',(I',)  VFs(x.r)  Then  Equation  (4)  becomes 


u  Vip(x,r) 


JO(x,t)|cos(y(x,t)-g>(x,f)) 


(5) 


We  model  the  target  as  a  collection  of  N  point  scatterers  and  write 

N 

P(y)~y£iA»e‘<>‘  5(y-y„|.  Each  of  the  independent  random  variables  {/!„}  and  {i>„} 

n=l 

are  assumed  to  be  uncorrelated  and  stnct-sense  stationary  for  all  n  =1,..  ,V  In  addition, 
we  assume  that  {  i>„ }  is  uniform  in  the  interval  (~n,n).  This  model  yields 


N=|X.S»I 


where  S„  h  Anexp{i(2kii ■  y„  +  i?„)] 


(6) 


(and  similarly  for  |G|)  It  is  easy  to  show  that,  for  large  N r  the  quantity  £  S„  is  a 
normally  distributed  mean  zero  random  variable  (this  follows  from  the  Central  Limit 

Theorem  [Reference  7]).  It  follows  that  |Fs(x,r)|  and  |G(x,r)|  are  independent  Rayleigh 


distributed  random  variables 


The  <p(x,t)  and  i/r(x,r)  are  independent  and  so  the  numerator  in  Equation  (5)  is  a 

normally  distributed  random  variable  The  probability  density  of  g(n)  =  u-V<p(x,r) 
follows  as  the  density  of  the  quotient  of  independent  normal  and  Rayleigh  random 
variables  We  obtain 
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(7) 


(8) 


The  target  shape  information  available  from  the  measured  {gj  comes  through  K2  which 
depends  upon  the  support  D  of  p(y)  and  can  be  deduced  by  comparing  the  statistics  of 
{gj  to  Equation  (7)  (see  below). 

The  simple  form  of  Equations  (7)  and  (8)  is  a  direct  consequence  of  our  target 

model — m  particular,  on  the  assumption  that  t?„  is  uniform  on  (-n,n).  However,  the 
qualitative  behavior  of  these  results  is  unaltered  by  more  general  assumptions  [Reference 
8).  We  have  opted  for  this  simpler  approach  because  it  is  significantly  more  tractable  and 
because,  at  high  frequencies  (<u  = 10  GHz),  the  net  error  is  equivalent  to  a  tiny  down- 
range  shift  in  p(y).  The  assumptions  on  D  are  more  important  since  if  the  support  of 

p(y)  is  too  restricted  then  the  sums  implied  by  the  integrals  (6)  will  not  satisfy  the 
requirements  of  the  central  limit  theorem.  For  this  reason,  we  consider  only  “complex” 
targets  (by  which  we  mean  targets  with  sufficient  scattering  centers  that  Equations  (8) 
and  (9)  are  good  approximations) 

One  additional  assumption  is  required.  In  practice,  the  direction  u  is  fixed  to  the 

receiver  and  so  the  statistics  of  {g,}  will  depend  on  target  orientation.  This  dependence 

is  not  nearly  so  sensitive  as  the  dependence  of  the  individual  measurements  g,  on  target 
aspect  (which  we  have  already  discussed)  However,  we  must  assume  that  the  target  is 

changing  its  aspect  at  a  sufficiently  slow  rate  that  the  data  set  {g,}  can  be  considered  to 
have  been  collected  over  a  small  domain  (say  <  10°)  of  target  orientations.  (In  practice, 
we  have  obtained  good  results  with  data  sets  of  only  a  few  dozen  samples  and  so,  this  is 
not  really  much  of  a  restriction.) 

The  derivative  of  the  scattered  field  phase  with  respect  to  wave  number  is  given  by 
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dtpjx.t) 


dk 


-  Im\ 


|e  i2 — 1 2"  yp(y)exp(,2kA  y)dy 

|£S(X.0|  D 


(9) 


Companng  (9)  with  Equation  (4),  we  can  see  that  the  statistics  of  this  measurement  are 
fundamentally  the  same  as  for  the  cross-range  derivatives  except  that  the  parameter 

ir2(n)  will  be  determined  by  the  down-range  properties  of  the  target  It  is  easy  to  see,  by 


substituting  u  =  i  cos  0  +  j  sin  6  into  Equation  (8)  that  k(9)  =  -\Jk2(0)  defines  an  ellipse 
whose  shape  and  orientation  depend  upon  the  cross-range  extent  of  the  target  The 

combination  ( ic2 (uj,  if2 (n))  will  define  an  ellipsoid  whose  shape  and  orientation  are 

determined  by  the  target  as  a  whole,  This  ellipsoid  can  be  mapped  to  overall  target 
dimensions  and  orientation.  However,  it  seems  doubtful  that  this  information  can  ever  be 
used  for  more  than  rough  target  categorization.  A  more  promising  approach  uses  a 

modified  estimate  k2 (u;t)  determined  as  a  function  of  the  target’s  down-range  extent 

The  wideband  case  The  narrowband  result  is  sufficiently  interesting  that  we  want  to 
preserve  as  much  of  its  character  as  possible.  To  do  this,  we  select  the  incident  signal  to 

be  of  the  form  F,(y,r)  =  £(n  y  -  ct)exp(ik0{n  y-ct))  where  k0  is  the  wave  number 

associated  with  a  “carrier”  frequency  co„  and  £(cr)  =  |E(cr)|e,ttn)  is  an  “envelope” 
function  with  compact  support. 

Let  £  =  2|x|  +  2n  y-ct.  We  require  £(g)  to  be  slowly  varying  in  comparison  with 
the  carrier  frequency: 


£'(£)| 


£(S) 


«K 


(10) 


Then  the  remaining  derivation  follows  as  before  and  we  can  conclude  that  /f(g,  u)  obeys 
Equation  (7)  with 
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**M  = 


(11) 


The  presence  of  |E(£)|2  in  (11)  reduces  the  support  of  the  summations  to  the 

intersection  of  the  support  of  \E(%)\  and  D  As  £  vanes  with  f,,  D  will  be  effectively 
“sliced”  into  a  collection  of  down-range  components.  The  statistics  of  (7)  and  (11)  allow 
us  to  associate  a  relevant  cross-range  target  structure  measure  to  each  down-range  slice 

ESTIMATION  OF  ; c2 

For  a  direction  u  fixed  to  the  receiver,  we  can  collect  a  set  of  phase-front  denvative 

data  {&},  l  =  1,.  .,P,  associated  with  P  (unknown)  target  aspect  angles.  Assuming  that 
these  aspects  angles  have  sufficiently  small  vanation  that  the  statistics  of  Equation  (7)  are 

quasi-stationary,  then  we  can  use  statistical  estimation  techniques  to  determine  ic2(u). 

The  maximum  likelihood  estimator  for  jc2(u)  will  be  the  one  that  maximizes  the 
function 


(12) 


(Here,  we  choose  to  maximize  the  log-likelihood  for  convenience.)  Differentiating 
In  (k-2)  with  respect  to  k2  and  setting  the  result  to  zero  yields 


£_y  J£VL=0 

1  *■«  „Z  ,  „2 

3  i=l  K  +S, 

Equation  (13)  is  the  sum  of  functions  of  the  form 


(13) 
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1  k>/2  (14) 

3  K2+g 2 


and  so  it  will  possess  ,,nly  two  real  solutions  which  are  located  symmetrically  about  the 
origin.  The  magnitude  of  either  zero  is  our  desired  estimate  and  can  be  obtained 
numerically  from  Equation  (13). 


THE  DOWN-RANGE  MAP  ASSOCLATED  WITH  A  t  Aivrs/i 

As  u  is  allowed  to  vary,  we  will  obtain  a  collection  of  estimate;,  tot  lf2(u,  —each  a 
solution  to  Equation  (13)  using  the  associated  input  data  {g,  =u  Vtpfx,,.);  In  a 
standard  Cartesian  coordinate  system,  with  axes  directions  i  and  j,  we  car.  associa  e  a 

polar  coordinate  angle  by  u  i  =  cos(6)  The  values  of  k(0)  will,  ideally,  trace  out  an 
ellipse  of  the  form 


K(0)  = 


d  _ 

l  +  ec<w[2(0  +  a)] 


(15) 


We  can  use  least-squares  estimation  to  fit  a  set  of  estimates  =  K2(0t)} ,  k  -  1,— £ 
to  the  ideal  form  of  Equation  (15)  This  will  allow  us  to  associate  the  parameter  set 

[d,  e,  a)  to  each  data  set  (V<p(x,,r)}.  These  parameters  are  related  to  target  cross-range 
extent  (through  <f),  cross-range  shape  (through  e),  and  target  orientation  (through  a) 
Moreover,  if  the  data  are  of  large  enough  bandwidth  that  we  can  isolate  separate  slices  of 
down-range  components  as  in  Equation  (11),  then  we  can  associate  these  structure  and 
onentauon  parameters  to  a  down-range  map  of  the  target. 


AN  EXAMPLE 


For  verification  purposes,  we  model  a  synthetic  target  as  a  collection  of  point 
scattered  with  fixed  amplitudes,  phases,  and  locations,  and  calculate  the  appropriate 
azimuth  and  elevation  data  that  would  be  obtained  by  a  phase-monopulse  tracking  system 
trained  on  the  randomly  maneuvering  target  We  have  created  these  data  with  sufficient 


baudwimii  diat  we  can  associate  i  =  1,  -P  data  values  to  each  of) 

down-range  target  slices.  These  data  are  vector-valued  with  components  consisting  of  the 
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azimuth  and  elevation  phase  derivatives.  The  interval  (0,  zr)  was  subdivided  into 
k  =  1,.  .,K  equal  subintervals  and  the  directions  u,  =  i  cos  9t  +  j  sin  ©,  were  used  to  create 
the  values  g,(t,)  =  ut  V<p(x,,f;)  which  were  input  into  Equation  (13)  to  estimate 
Finally,  these  values  were  fit  to  the  ellipses  of  Equation  (15)  to  determine  the 
target  parameters  {d(r;  j,  e(t, ),  «(f, )}  • 

Figure  1  displays  a  collection  of  200  point  scatterers  randomly  placed  on  the  skin  of  a 
synthetic  target.  From  this  collection  of  scatterers  we  calculated  the  (wideband)  azimuth 
and  elevation  phase  monopulse  tracking  data  for  each  of  100  random  target  orientations 
cnosen  from  within  a  10°  x  10°  domain  centered  on  the  target’s  axis.  We  have  modeled 
the  incident  signal  as  positively  directed  along  the  z-axis  with  carrier  frequency  30  GHz 
and  simple  unit-box  envelope  characterized  by  a  bandwidth  of  1  2  GHz.  Figure  2  shows 
the  results  of  the  algorithm  displayed  in  graphical  format  This  image  is  the  “enhanced 
down-range  map”  and  clearly  shows  the  influence  of  target  cross-range  structure  in  an 
interpretable  way. 


FIGURE  1.  Sample  Target  Plotted  are  the  location  of  the  200  scattering  centers 
in  space  These  scattenng  centers  were  chosen  to  be  randomly  located  on  the 
surface  of  the  target  and  obey  the  approximate  relation  of  10  scatterers  per  down- 
range  slice  (Target  support  and  onentation  are  shown  m  the  inset )  Not  plotted 
are  the  intensity  and  phase  shift  associated  with  each  scatterer. 
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(a)  Displa;  s  the  down-range  map  as  estimated  by  the  algorithm 


FIGURE  2  Enhanced  Down-Range  Maps  Associated  With  Target  1. 


Note  that  errors  in  a  diminish  the  interpretability  of  the  visual  reconstruction.  An 
easy  enhancement  which  helps  to  resolve  this  ellipse  misalignment  is  to  average  a 
overall  down-range  slices  for  which  the  ratio  of  semi-major  to  semi-minor  axes  exceeds 
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seme  threshold  value.  We  have  not  performed  this  averaging  here  in  order  that  only  the 
relevant  properties  of  the  basic  algorithm  be  displayed  The  down-range  reconstruction 
of  Figure  2  required  1  2  GHz  of  bandwidth  to  achieve  the  displayed  resolution  Less 
bandwidth  results  in  less  resolution — 600  MHz  of  bandwidth,  for  example,  would  allow 
us  to  plot  only  half  of  the  ellipses  of  Figure  2.  However,  in  our  studies  we  have  seen  that 
even  bandwidths  as  small  as  100  MHz  have  resulted  in  meaningful  classification  of  10 
meter  (synthetic)  targets  while  standaid  (simple)  down-range  profiles  have  proven  to  be 
essentially  useless  Even  for  narrowband  systems  (vanishing  bandwidth)  we  can  still 
determine  target  size  and  orientation 

The  example  target  not  is  intended  to  represent  the  actual  distribution  of  scattering 
centers  for  any  particular  aircraft  Real  targets  cannot  be  expected  to  have  a  uniform 
distribution  of  scatterers  lying  on  the  surface  of  a  given  cross-sectional  slice  and  while 
they  can  be  quite  complex,  they  may  not  be  composed  of  as  many  as  200  points.  (Or 
models  were  simple  expedients  and  the  total  number  of  scattering  centers  was  chosen  by 
requiring  =10  scatterers  per  slice )  However,  a  particular  type  of  aircraft  can  be  expected 
to  possess  a  unique  distribution  of  scatterers  and  by  virtue  of  the  squared  radius  weight 
they  are  given  in  Equation  (11),  the  current  algorithm  should  prove  to  be  a  sensitive 
classifier.  Furthermore,  the  complexity  of  the  aircraft  (i  e.,  the  number  of  scattering 
centers)  will  generally  increase  with  increasing  earner  frequency  co0  and  the  number  of 
scatters  per  slice  can  approach  the  example  values  under  a  wide  vanety  of  situations.  (A 
detailed  examination  of  the  dependence  of  algonthm  accuracy  on  scatterer  number  is 
beyond  the  scope  of  the  present  work  Some  initial  work  along  these  lines  was  attempted 
in  Reference  5 ) 


DISCUSSION  AND  CONCLUSION 


The  phase-front  of  the  field  scattered  from  a  complex  target  contains  a  great  deal  of 
information  about  that  target  In  the  forgoing  discussion,  we  have  shown  how  some  of 
that  information  can  be  retrieved  by  ordinary  missile  seeker  systems.  In  addition,  we 
have  detailed  a  method  that  relates  cross-range  target  structure  to  a  set  of  parameters 
which  can  be  incorporated  into  a  machine-based  target  classification  system. 

The  proposed  technique  has  been  tested  using  both  real  and  synthetic  data  and  has 
proven  to  be  an  effective  discriminator  of  overall  target  size  and  general  shape.  In 
addition,  the  method  has  continued  to  yield  useful  informauon  about  the  target  even  when 
the  data  are  few  in  number  and  the  model  assumptions  are  significantly  relaxed  Of 
course,  throughout  this  discussion  we  have  implicitly  assumed  that  the  radar  is  “trained 
on  the  target.”  Ir.  practice,  however,  the  so-called  '‘glint  problem”  may  prevent  this  from 
being  done  with  complete  accuracy 

The  glint  problem  (also  known  as  angular  scintillation!  refers  to  the  error  between  the 
target  bearing  predicted  by  the  direction  of  the  phase  gradient  and  the  true  target  beanng. 
This  difference  is  also  the  data  that  we  are  using  for  our  cross-range  structure  algonthm 
and  cannot  be  estimated  if  the  true  target  beanng  is  unknown.  Significant  progress  has 
been  made  in  solving  this  problem  and  in  fact,  it 's  known  that  tracking  error  can  be 
effectively  eliminated  by  correctly  averaging  the  separate  frequency  domain  returns 
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collected  by  wideband  systems  [Reference  9]  It  is  precisely  these  wideband  radar 
systems  that  we  have  tailored  this  algorithm  to.  (Some  progress  has  also  been  made  in 
eliminating  tracking  error  from  narrowband  systems  [Reference  10]  and  the  ordinary 
narrowband  version  of  the  algorithm  can  be  expected  to  be  implernentable  as  well ) 


Future  plans  include  a  more  complete  verification  of  the  method  using  a  larger  variety 
of  measured  data  In  addition,  we  hope  to  investigate  the  classifiability  of  specific 
aircraft  targets  by  their  phase-front  derivatives 
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